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Abstract 

Within the four-loop e expansion, we study the critical behavior of certain 
antiferromagnets with complicated ordering. We show that an anisotropic stable 
fixed point governs the phase transitions with new critical exponents. This is 
supported by the estimate of critical dimensionality A'^^ = 1.445(20) obtained 
from six loops via the exact relation = established for the real and 

complex hypercubic models. 
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It is known that the critical properties of phase transitions in certain antiferromag- 
nets involving an increase of the unit cell in one or more directions at the critical tem- 
perature can be described by a generalized 2A^-component {N > 2) Ginzburg-Landau 
model with three independent quartic terms 
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associated with the isotropic, cubic, and tetragonal interactions, respectively |1| . Here 
ifi is the real vector order parameter in D = 4 — e dimensions and itlq is proportional 
to the deviation from the mean-field transition point. When N = 2, Hamiltonian 
(|l]) describes the antiferromagnetic phase transitions in TbAu2 and DyC2 and the 
structural phase transition in Nb02 crystal^. Another physically important case = 3 
is relevant to the antiferromagnetic phase transitions in such substances as K2lrCl6, 
TbD2, MnS2, and Nd. All these phase transitions are known from experiments to 
be of second orderQ (see Ref. and references therein). However, the experimental 
data were insufficiently accurate to provide reliable values of critical exponents and 
the obtained estimates |0, ||, ^ were found to differ significantly from the theoretically 
expected numbers. 

For the first time the magnetic and structural phase transitions described by model 
(|l]) were studied in the framework of the renormalization group (RG) by Mukamel and 
Krinsky within the lowest orders in e: ^. A three-dimensionally stable fixed point 
(FP) with coordinates u* > 0, v* = z* > was predicted^ That point was shown to 
determine a new universality class with a specific set of critical exponents. However, 
for the physically important case N = 2, the critical exponents of this unique stable 
FP turned out to be exactly the same as those of the 0(4)symmetric one. 

For the years an alternative analysis of critical behavior of the model, the RG 
approach in three dimensions, was carried out within the two- and three-loop approxi- 



mations [0, |Tl|. Those investigations gave the same qualitative predictions: the unique 
stable FP does exist on the 3D RG fiow diagram. By using different resummation pro- 
cedures, the critical exponents computed at this point proved to be close to those of 
the Bose FP [u = 0, v = z > 0) rather than the isotropic (O(A^)-symmetric; u > 0, 
V = z = 0) one. It was also shown that the unique and the Bose FPs are very close 
to each other, so that they may interchange their stability in the next orders of RG 
approximation |TT]. 



""^The phase transitions in helical magnets Tb, Dy, and Ho belong to the same class of universality 

^An interesting type of multisublattice antiferromagnets, such as MnO, CoO, FeO, and NiO, was 
studied in Ref. |, |, I- It was shown, in the leading orders in e, that the phase transitions in these 
substances are of first order. 

•^Following Mukamel 0, we call this point "unique". 
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Recently, the critical properties of the model were analyzed in third order in e [|T2 



TB|. Investigation of the FP stability and calculation of the critical dimensionality A''^ 
of the order parameter, separating two different regimes of critical behavioif], confirmed 
that model ([I|) has the unique stable FP at = 2 and = 3. However, the twofold 
degeneracy of the stability matrix eigenvalues at the one-loop level was observed for this 
FP |T^. That degeneracy was shown to cause a substantial decrease of the accuracy 



expected within the three-loop approximation and powers of yfe to appear in the 
expansions^. So, computational difficulties were shown to grow faster than the amount 
of essential information one may extract from high-loop approximations. That resulted 
in the conclusion that the 5-expansion method is not quite effective for the given model. 

Another problem associated with model (P is the question whether the unique FP 
is really stable in 3D, thus leading to a new class of universality, or its stability is only 
an effect of insufficient accuracy of the RG approximations used. Indeed, there are 
general nonperturbative theoretical arguments indicating that the only stable FP in 
3D may be the Bose one and the phase transitions of interest should be governed by 
that stable FP [1^. However, up to now this assertion found no confirmation within 



the RG approach. In such a situation it is highly desirable to extend already known 
e expansions for the stability matrix eigenvalues, critical exponents, and the critical 
dimensionality in order to apply more sophisticated resummation technique to longer 
expansions. 

In this Letter we, firstly, avoid the problem of the eigenvalues degeneracy in model 
(|1]) by analyzing the critical behavior of an equivalent complex A^*-^-component order 
parameter model with the effective Hamiltonian 



comprising the isotropic and cubic interaction^. Note that this Hamiltonian also 
describes the real hypercubic model ll^ if ijji is thought to be the real A^^-component 



order parameter. The model (|^) comes out exactly from model (|T]) at vq = zq and 
it is free from the eigenvalues degeneracy. Secondly, we examine the existence of the 
anisotropic stable FP in model (|^) on the basis of the higher-order e expansion. Namely, 
using dimensional regularization and the minimal subtraction scheme |]18[, we derive 



the four-loop RG functions as power series in e and analyze the FP stability. For 
the first time, we give realistic numerical estimates for the stability matrix eigenvalues 
using the Borel transformation with a conformal mapping [0. This allows us to carry 
out the careful analysis of the stability of all the FPs of the model. We state the exact 
relation A"^ = between the critical (marginal) spin dimensionalities of the real 
and complex hypercubic models and obtain the estimate A^^ = 1.445(20) using six- 
loop results of Ref. |2^. We show that the anisotropic (complex cubic; u 0, v ^ 0) 



"^When N > Nc the unique FP is stable in 3D while for N < Nc the stable FP is the isotropic one. 

^Similar phenomenon was observed earlier in studying the impure Ising model (see Refs. fl^ ). 
Half-integer powers in s arising in that model have different origin but also lead to the loss of accuracy. 

^The model with the complex vector order parameter was considered by Dzyaloshinskii |l6|| in 
studying the phase transitions in DyC2, TbAu2 (N^ = 2) and TbD2, MnS2, and Nd {N^ = 3). 
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stable FP of model @, being the counterpart of the unique point in model (|T]), does 
exist on 3D RG flow diagram at A^"-^ > A*"^. For this stable FP we give more accurate 
critical exponents estimates in comparison with the previous three-loop results [0 by 
applying the summation technique of Ref. to the longer series. 

The four-loop e expansion for the /3-functions of model (|^) were recently obtained 
by the present authors in Ref. [Q. From the system of equations f3u{u*,v*) = 0, 
(3v{u*, V*) = one can calculate formal series for the four FPs: the trivial Gaussian one 
and nontrivial isotropic, Bose, and complex cubic FPs. Instead of presenting here the 
FPs themselves, which have no direct physical meaning, we present the eigenvalues of 
the stability matrix 

C dt3u{u,v) df3u{u,v) \ 
df3%,v) d(3%,v) j (3) 

taken at the most intriguing Bose and complex cubic FPs. They are 



1 6 2 1 

2 20 8 



257 384 



1 

+ I6 



5152 + 5!ic(3) 

1250 125^ ' 



125 125 
576 



C(3) 



125 



C(4,.1^C.) 



^2 



1 


14 




10 


100 


1 


61 






+ 16 


250 



5^ + 



311 768 



625 625 



C(3) 



3125 

at the Bose FP and for = 2 

1 13 2 
-1 = -^e^-e + 



625 



625 



65 
36 



rC(3) 



1 

+ 16 



(4) 



and A^^ = 3 



CU2 



12 8 



- + -C(3) 
18 6^^ ^ 



+ 



16 



181 



50 



145 5 
^-^C(3) + -C(4)-7C(5) 



(5) 



1 58 2 1 
= "2^+ 220^ +8 



+ 



19533 14832 



1 

16 



310518757 



15125 
1644864 



C(3) 



91506250 1830125 



C(3) 



6655 
22248 



6655 



283056^,^; 



3 



0^2 



1 2 

~22^ ^ 2420 



90363 3408 



1 

+ 16 



1151231173 



166375 
50696504 



C(3) 



1006568750 20131375 



C(3) 



73205 
5112 



73205 



107344^,^, 
161051^^^) 



(6) 



at the complex cubic one, where C(3)) C(4), and C(5) are the Riemann C, functions. 

It is known that RG series are at best asymptotic. An appropriate resummation 
procedure has to be apphed in order to extract rehable physical information from 
them. To obtain the eigenvalue estimates we have used an approach based on the 
Borel transformation modified with a conformal mapping |jl9|, ^ . If both eigenvalues 
of matrix (0) are negative, the associated FP is infrared stable and the critical behavior 
of experimental systems undergoing second-order phase transitions is determined only 
by that stable point. For the Bose and the complex cubic FPs our numerical results are 
presented in Table I. It is seen that the complex cubic FP is absolutely stable in D = 3 
{e = 1), while the Bose point appears to be of the "saddle" type. However u!2's of 
either points are very small at the four-loop level, thus implying that these points may 
swap their stability in the next order of RG approximation. We can compare uj2 at the 
complex cubic FP quoted in Table I with the three-loop results of Ref. 
in the framework of RG approach directly in 3D. Those estimates U2 
= 2 and uj2 = —0.011 for A^*^ = 3 are solidly consistent with ours. 

The four-loop e expansion for the critical dimensionality of the order parameter of 
model f0) reads 



10 



obtained 
0.010 for 



n: 



c 



2-e + 



24 



6C(3) - 1 



1 



144 



45C(3) + 135C(4) - 600C(5) - 1 



Instead of processing this expression numerically, we state the exact relation = 
\N^, which is independent on the order of approximation used. In fact, the critical 
dimensionality for the complex cubic model is determined as that value of iV", at 
which the complex cubic FP coincides with the isotropic one. In the same way, the 
critical dimensionality A*"^ is defined for the real cubic model. Both systems exhibit 
effectively the isotropic critical behavior at A^*" = and A^'^ = N^. So, because 
the complex 0(2A^'^)-symmetric model is equivalent to the real 0(A^^) -symmetric one, 
the relation 2N'^ = holds true. For A^*^ > the complex cubic FP of model (^) 
should be stable in 3D. 

The five-loop e expansion for A^ was recently obtained in Ref. 
tion of that series gave the estimate A^f = 2.894(40) (see Ref. p4 
conclude that A^ = 1.447(20) from the five-loops. Practically the same estimate 
A^ = 1.435(25) follows from a constrained analysis of A^ taking into account A*"^ = 2 
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Resumma- 
Therefore we 



From the recent pseudo-e expansion analysis of the real hyper- 
one can extract A5 = 1.431(3). However the most accurate estimate 



in two dimensions PO 
cubic model 

= 1.445(20) results from the value A'^^ = 2.89(4) obtained on the basis of the 
numerical analysis of the four-loop and the six-loop 3D RG expansions for the 
/3-functions of the real hypercubic model. 
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Finally, we have computed the four-loop e series for the critical exponents. At the 
stable complex cubic FP they are 
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for A^*^ = 3. Other critical exponents can be found through the known scaling relations. 
The numerical estimates obtained are collected in Table II. The critical exponents for 
the isotropic and the Bose FPs are also presented, for comparison. We can compare our 
results with the available experimental data. For example, in the case of the structural 
transition in the Nb02 crystal the critical exponent of spontaneous polarization was 
measured in Ref. P], 0.33 < P < 0.44. Our estimate /3 = 0.371 obtained using the 
data of Table II and scaling relations lies in that interval. 

In summary, the four-loop e-expansion analysis of the GinzburgLandau model with 
cubic anisotropy and complex vector order parameter relevant to the phase transitions 
in certain antiferromagnets with complicated ordering has been carried out. Investi- 
gation of the global structure of RG flows for the physically significant cases N'" = 2 
and A^*" = 3 leads to the conclusion that the anisotropic complex cubic FP is abso- 
lutely stable in 3D. Therefore the critical thermodynamics of the phase transitions in 
the Nb02 crystal and in the antiferromagnets TbAu2, DyC2, K2 IrClg, TbD2, MnS2, 
and Nd should govern by this stable point with a specific set of critical exponents, in 
the frame of the given approximation. The critical dimensionality = 1.445(20) 
obtained from six loops supports this conclusion. At the complex cubic FP, the crit- 
ical exponents were calculated using the Borel summation technique in combination 
with a conformal mapping. For the structural phase transition in Nb02 and for the 
antiferromagnetic phase transitions in TbAu2 and DyC2 they were shown to be close 
to the critical exponents of the 0(4)-symmetric model. In contrast to this, the critical 
exponents for the antiferromagnetic phase transitions in K2lrCl6, TbD2, MnS2, and Nd 
turned out to be close to the Bose ones. 
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Table I: Eigenvalue exponent estimates for the Bose (BFP) and the complex cubic 
(CCFP) FPs at A^*^ = 2 and A^*" = 3 obtained in the four-loop order in e {e = 1) using 
the Borel transformation with a conformal mapping. 



Type of 
FP 
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002 
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-0.400(30) 
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Table II: Critical exponents for the isotropic (IFP), the Bose (BFP), and the complex 
cubic (CCFP) FPs at A^*" = 2 and A^*" = 3 calculated in the four-loop order in e 
{e = 1) using the Borel transformation with a conformal mapping. 



Type of 
FP 


A^^' = 2 


N^ = 3 


V 


u 


7 


V 


u 
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IFP 


0.0343(20) 


0.725(15) 


1.429(20) 


0.0317(10) 


0.775(15) 


1.524(25) 


BFP 


0.0348(10) 


0.664(7) 


1.309(10) 


0.0348(10) 


0.664(7) 


1.309(10) 


CCFP 


0.0343(20) 


0.715(10) 


1.404(25) 


0.0345(15) 


0.702(10) 


1.390(25) 



Although our calculations show that the complex cubic FP, rather than the Bose 
one, is stable at the four-loop level, the eigenvalues 002 of both points are very small. 
Therefore the situation is very close to marginal, and the FPs might change their stabil- 
ity to opposite in the next order of perturbation theory, so that the Bose point would 
occur stable. This conjecture is in agreement with the recent six-loop RG study of 
three-coupling-constant model (|l|) directly in three dimensions p6| . The authors argue 
the global stability of the Bose FP, although the numerical estimate 002 = —0.007(8) 
of the smallest stability matrix eigenvalue of the Bose point appears to be very small 
and the apparent accuracy of the analysis does not exclude the opposite sign for 002. 
In this situation it would be very desirable to compare the critical exponents values 
obtained theoretically with values that could be determined from experiments, in order 
to verdict which of the two FPs is really stable in physical space. Finally, it would be 
also useful to investigate certain universal amplitude ratios of the model because they 
vary much more among different universality classes than exponents do and might be 
more effective as a diagnostic tool. 

We are grateful to Professor M. Henkel for helpful remarks and to Dr. E. Vicari 
for sending to us a copy of the review cited in Ref. [^. This work was supported 
by the Russian Foundation for Basic Research via grant No. 01-02-17048, and by the 
Ministry of Education of Russian Federation via grant No. EOO-3. 2-132. 



6 



References 

[1] D. Mukamel, Phys. Rev. Lett. 34, 481 (1975); D. Mukamel and S. Krinsky, J. 
Phys. C 8, L496 (1975); Phys. Rev. B 13, 5065 (1976). 

[2] P. Bak and D. Mukamel, Phys. Rev. B 13, 5086 (1976). 

[3] S. A. Brazovskii and I. E. Dzyaloshinskii, Pis'ma Zh. Eksp. Teor. Fiz. 21, 360 

(1975) [JETP Lett. 21, 164 (1975)]; S. A. Brazovskii, I. E. Dzyaloshinskii, and B. 
G. Kukharenko, Zh. Eksp. Teor. Fiz. 70, 2257 (1976) [Sov. Phys. -JETP 43, 1178 

(1976) ] and references therein. 

[4] P. Bak, S. Krinsky, and D. Mukamel, Phys. Rev. Lett. 36, 52 (1976). 
[5] D. Mukamel and S. Krinsky, Phys. Rev. B 13, 5078 (1976). 

[6] J. C. Toledano, L. Michel, P. Toledano, and E. Brezin, Phys. Rev. B 31, 7171 
(1985). 

[7] J. Als-Nielsen, in: Phase transitions and Critical Phenomena, edited by C. Domb 
and M. S. Green (Academic, New York, 1976), Vol. 5A, p. 88. 

[8] R. Pynn and J. D. Axe, J. Phys. C 9, L199 (1976). 

[9] J. Schneck et al., Phys. Rev. B 25, 1766 (1982). 

[10] N. A. Shpot, Phys. Lett. 133A, 125 (1988); ibid 142A, 474 (1989). 

[11] K. B. Varnashev and A. I. Sokolov, Fiz. Tverd. Tela (St. Petersburg) 38, 3665 
(1996) [Phys. Sohd State 38, 1996 (1996)]; A. I. Sokolov, K. B. Varnashev, and 
A. I. Mudrov, Int. J. Mod. Phys. B 12, 1365 (1998); A. I. Sokolov and K. B. 
Varnashev, Phys. Rev. B 59, 8363 (1999). 

[12] K. De'Bell and D. J. W. Geldart, Phys. Rev. B 32, 4763 (1985); E. J. Blagoeva, 
Mod. Phys. Lett. B 10, 439 (1996). 

[13] A. I. Mudrov and K. B. Varnashev, Phys. Rev. B 57, 3562; ibid 5704: (1998). 

[14] D. E. Khmclnitskii, Zh. Eksp. Teor. Fiz. 68, 1960 (1975) [Sov. Phys. JETP 41, 
981 (1976)]. 

[15] J. Sak, Phys. Rev. B 10, 3957 (1974); K. A. Cowley and A. D. Bruce, J. Phys. C 
11, 3577 (1978). 

[16] I. E. Dzyaloshinskii, Zh. Eksp. Teor. Fiz. 72, 1930 (1977) [Sov. Phys. -JETP 45, 
1014 (1977)]. 

[17] A. Aharony, in Phase Transition and Critical Phenomena, edited by C. Domb and 
M. S. Green (Academic Press, New York, 1976), Vol. 6, p. 357. 



7 



[18] G. 't Hooft and M. Veltman, Nucl. Phys. B 44, 189 (1972); G. 't Hooft, Nucl. 
Phys. B 61, 455 (1973). 

[19] J. C. Le Guillou and J. Zinn- Justin, Phys. Rev. Lett. 39, 95 (1977); Phys. Rev. B 
21, 3976 (1980); R. Seznec and J. Zinn- Justin, J. Math. Phys. 20, 1398 (1979); A. 
A. Vladimirov, D. I. Kazakov, and O. V. Tarasov, Zh. Eksp. Teor. Fiz. 77, 1035 
(1979) [Sov. Phys. -JETP 50, 521 (1979)]. 

[20] J. M. Carmona, A. Pehssetto, and E. Vicari, Phys. Rev. B 61, 15136 (2000). 

[21] A. I. Mudrov and K. B. Varnashev, Phys. Rev. E 58, 5371 (1998). See also in: 
Problems of Quantum Field Theory (Dubna: JINR, 1999), p. 267; preprint hep- 
I th/981lT25| . 

[22] A. I. Mudrov and K. B. Varnashev, J. Phys. A 34, L347 (2001). See also in: Group 
Theoretical Methods in Physics (Dubna: JINR, 2001); preprint |cond- mat / 1 1 1 6^ . 

[23] H. Kleinert, V. Schulte-Frohlinde, Phys. Lett. 342B, 284 (1995). 

[24] K. B. Varnashev, J. Phys. A 33, 3121 (2000); Phys. Rev. B 61, 14660 (2000). 

[25] R. Folk, Yu. Holovatch, and T. Yavors'kii, Phys. Rev. B 62, 12195 (2000). 

A. Pelissetto and E. Vicari, preprint |cond-mat / 1 2 1 64 . 



8 



